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Abstract Many inversion-layer metal-insulator-semiconductor (IL/MIS) solar cells have been 
fabricated. There are around eighteen 1 cm 2 IL/MIS solar cells which have efficiencies greater 
than 7%. There are only about three 19 cm 2 IL/MIS cells which have efficiencies greater than 
4%. The more accurate control of the thickness of the thin layer of oxide between aluminum and 
silicon of the MIS contacts has been achieved. A lot of effort and progress have been made in this 
area. 

A comprehensive model for MIS contacts under dark conditions has been developed that 
covers a wide range of parameters. It has been applied to MIS solar cells. One of the main 
advantages of these models is the prediction of the range of the thin oxide thickness versus the 
maximum efficiencies of the MIS solar cells. This is particularly important when the thickness is 
increased to 25A. This study is very useful for our investigation of the EL/MIS solar cells. 

The two-dimensional numerical model for the IL/MIS solar cells has been tried to develop 
and the results are presented in this report. 

1.0 Introduction 

The objectives of this study are: (1) investigation and fabrication of large area silicon 
inversion-layer metal-insulator-semiconductor (MIS) solar cells for use in space solar arrays, (2) 
continue developing a simple inexpensive process for fabricating high-efficiency IL solar cells, and 
(3) continue modeling the IL/MIS solar cells. 

Man y inversion layer MIS solar cells have been fabricated. There are around eighteen 
1 cm 2 IL/MIS solar cells which have efficiencies greater than 7%. There are only about three 

19 cm 2 IL/MIS cells which have efficiencies greater than 4%. 

The growth of the thin layer of oxide between aluminum and silicon of the MIS contacts 
was one of the major tasks in this project. A lot of effort and progress have been made in this 
area [1,2]. The theoretical study of MIS contacts and MIS solar cells have been pursued. A 
comprehensive model for MIS devices has been developed. It has been applied to MIS solar cells 
[3,4], This study is very useful for our investigation of the inversion-layer MIS solar cells. 
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The two-dimensional numerical model for the IL/MIS solar cells has been tried to develop, 
but it fails to yield good results. 

2.0 The Fabrication of IL/MIS Solar Cells 

The operational principle, the design consideration, and the fabrication of inversion-layer 
MIS solar cells have been presented in other papers [5, 6, 7,8]. The fabrication procedures of the 
inversion-layer MIS solar cells are briefly described as follows: 

1. Silicon dioxide (Si0 2 ) with a thickness of approximately 1300A is thermally grown on a 
<1 1 1> p-type silicon substrate. 

2 The front metal contact regions are defined using the first mask. They are etched on the front 
surface. 

3. A thin oxide layer is grown in metal contact regions. 

4. Aluminum is evaporated on the back surface of the silicon substrate and followed by sintering 
in a 470°C furnace. 

5. Aluminum is evaporated on the front side. 

6. The grid patterns of the front totalization are defined by using a second mask. 

We have fabricated many IL/MIS solar cells. The 1 cm 2 cells with 7% or higher 
efficiencies are listed in Table 1. The 1 cm 2 cells with 6% to 7% efficiencies are tabulated in 
Table 2. The 1 cm 2 IL/MIS solar cells with efficiencies between 4% and 6% are shown in Table 
3. The large area (total area — 19.78 cm 2 ) IL/MIS solar cells are summarized in Table 4 [2]. 
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Table 1. 1 CM 2 IL-MIS Solar Cells with 7% or Higher Efficiencies 


Substrate 

Resistivity 

Q-cm 

Cell ID 

Experimental Results 

V 0C 

(V) 

^sc 

(mA) 

FF 

% EFF 

1.2 

8-4 

H 

0.567 

29.050 

0.444 

7.003 

1.2 

1-1 

H 

0.575 

30.400 

0.536 

8.909 

B 

0.570 

29.000 

0.452 

7.106 

1.2 

1-1 

H 

0.536 

27.350 

0.561 

7.825 


0.526 

27.100 

0.552 

7.481 

0.7-1. 3 

13-1 

H 

0.562 

26.200 

0.546 

7.646 

B 

0.559 

26.550 

0.552 

7.787 

13-2 

H 

0.555 

26.600 

0.560 

7.863 

B 

0.556 

25.900 

0.548 

7.500 

C 

0.548 

25.300 

0.561 

7.402 

D 

0.560 

27.050 

0.582 

8.383 

0.7-1. 3 

14-2 

H 

0.549 

24.600 

0.578 

7.430 

3. 0-6.0 

14-3 

H 

0.525 

25.300 

0.589 

7.441 

B 

0.522 

24.750 

0.577 

7.088 

C 

0.526 

25.100 

0.590 

7.404 

0.7-1. 3 

2-2 

C 

0.547 

26.400 

0.614 

8.430 

0.7-1. 3 

4-1 

H 

0.523 

25.300 

0.564 

7.100 

mm 

c 

0.536 

25.850 

0.612 

8.067 
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Table 2. 1 CM 2 EL-MIS Solar Cells with 6% to 7% Efficiencies 


Experimental Results 


Substrate 

Resistivity 

Q-cm 

Cell ID 

V 0C 

(V) 

Isc 

(mA) 

FF 

% EFF 


8.3 

m 

0.573 

29.000 

0.406 

6.456 

1.2 


B 

0.566 

27.700 

0.432 

6.477 


8.4 

C 

0.557 

27.450 

0.453 

6.628 

1.2 

9.3 

m 

0.579 

30.350 

0.376 

6.331 



c 

0.575 

29.700 

0.387 

6.319 

0.7-1. 3 

13-1 

c 

0.563 

25.600 

0.459 

6.291 



D 

0.563 

25.550 

0.464 

6.350 



B 

0.546 

24.750 

0.544 

6.994 

0.7-1. 3 

14-2 

C 

0.537 

24.250 

0.530 

6.563 



D 

0.545 

24.550 

0.539 

6.865 

3. 0-6.0 

14-3 

D 

0.524 

23.800 

0.577 

6.849 



B 

0.502 

26.750 

0.496 

6.333 

8.0-12.0 

14-5 

D 

0.507 

26.900 

0.506 

6.565 



H 

0.494 

22.650 

0.608 

6.476 

3. 0-6.0 

1-3 

D 

0.496 

22.950 

0.645 

6.986 


2-1 

R9 

0.523 

24.800 

0.546 

6.737 

0.7-1. 3 


c 

0.525 

23.950 

0.524 

6.263 


2-2 

D 

0.534 

23.600 

0.543 

6.511 










































































Table 3. 1 CM 2 II^MIS Solar Cells with 4% to 6% Efficiencies 


Substrate 

Resistivity 

Q-cm 

Cell ID 

Experimental Results 

V 0C 

(V) 

Isc 

(mA) 

FF 

% EFF 

1.2 

9-1 

mm 

0.566 

29.200 

0.317 

5.008 

9-2 

mm 

0.557 

28.050 

0.353 

5.275 

c 

0.555 

27.450 

0.324 

4.730 

0.7-1. 3 

1-1 

B 

0.491 

24.400 

0.475 

5.412 

c 

0.482 

20.700 

0.451 

4.280 

0.7-1. 3 

3-1 

■■ 

0.536 

26.450 

0.421 

5.674 

B 

0.525 

25.000 

0.417 

5.212 

3-2 


0.531 

25.550 

0.407 

5.255 

c 

0.537 

25.400 

0.354 

4.592 

0.7-1. 3 

mm 

D 

0.542 

28.300 

0.330 

4.811 

8.0-12.0 

2-3 

■« 

0.496 

27.050 

0.404 

5.150 

K 

0.505 

28.100 

0.331 

4.463 

0.7-1. 3 

92-2 

nn 

0.533 

27.750 

0.357 

4.483 

0.7-1. 3 

2 

c 

0.507 

23.350 

0.418 

4.704 

D 

0.505 

24.100 

0.409 

4.734 

0.7-1. 3 

4-1 

C 

i 

0.524 

23.050 

0.379 

4.350 

mm 

B 

0.535 

24.500 

0.402 

5.008 


5 











































































Table 4. 19 CM 2 IL-MIS Solar Cells 


Substrate 

Resistivity 

O-cm 

Cell ID 

Thin 

Oxide 

Thickness 

(Average) 

Experimental Results 

V 0C 

(V) 

Isc 

(mA) 

FF 

% EFF 

0.7-1. 3 

1-H 

22k 

0.524 

342.5 

0.305 

2.438 

1-WO 

19A 

0.481 

171.0 

0.473 

1.733 

1-M 

19A 

0.452 

128.6 

0.466 

1.208 


2-H 

2lA 

0.524 

210.0 

0.401 

1.968 

0.2-0. 5 

2-WO 

19A 

0.482 

153.6 

0.447 

1.476 


2-M 

18A 

0.404 

70.8 

0.512 

0.653 

3. 0-6.0 

; 3-H 

2lA 

0.540 

517.0 

0.323 

4.014 

3-WO 

19A 

0.536 

571.0 

0.394 

5.374 

3-M 

18A 

0.531 

534.0 

0.384 

4.855 


3.0 A Comprehensive Analytical Model for Metal-Insulator-Semiconductor 
(MIS) devices 

A detailed model for MIS devices has been developed that covers a wide range of 
parameters which include surface states, silicon dioxide thickness, substrate doping, fixed oxide 
charges, substrate thickness, and metal work function [3,4]. One of the main advantages of this 
model is the prediction of the range of the thin oxide thickness versus the maximum efficiencies of 
the MIS solar cells. In this model we have made the following three main assumptions: 

1 . The semiconductor charges Q sc is calculated under non-equilibrium conditions, which can be 

written as: 
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Qsc = 


P'P, +e~^ -l + -^-{^[^-l] -P^} 


( 1 ) 


The old expression to calculate Q sc [9] assumes equilibrium condition for Q sc It gives 
expression like equation (1), but without the term e^ J . 


Qsc = +JlKT^N~ A 


P %+e 






( 2 ) 


2. The tunneling factor ( 0(x,8) = exp{J \K\dx ) ) is calculated from the curves generated by Card 
[10,1 1] for J|AT| dx = y} 12 ^ against . 

3 . The permittivity of oxide is taken to be that of free space. 

Without the first and second assumptions we could predict that maximum efficiencies 
0 0 0 0 
occur at about 10 A <d t ■ < 15 A (the range of the thin oxide thickness is from 15 A to 25 A) 

and that the efficiency is very sensitive to the variation of the oxide thickness. If the efficiency is 

so sensitive as the old model claimed, higher efficiencies could not be obtained in practice. Our 

model includes the first and second assumption and predicts more practical results, which are 

O 

shown in Table 5 and Figures 1 and 2. From Table 5 we see that for d i ■ = 15 A, the efficiencies 
with and without the first and second assumptions are nearly the same. That is why the old 
models can predict the efficiency correctly at this oxide thickness. When the thickness of the thin 

O 

oxide is increased to 20A our model gives 17.3% efficiency instead of 14.1% calculated without 

O 

the first and second assumptions. If the thickness is increased to 25 A our model predicts 17.1% 
efficiency which is compared with 8.76% calculated without these assumptions. In our model the 

o 

efficiency starts decreasing only with d t =25 A. This is clearly indicated in Figure 1. Figure 3 
shows the 1-V characteristic of MIS solar cell for different oxide thickness with and without 
these assumptions. For the thin oxide in MIS contacts our model for MIS devices predicts more 
practical results [10,12,13,14], This is particularly important when the thickness is increased to 
25 A [4,14], 
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Table 5. Effects of the New Assumptions in our Model 


Eq. # 

d; (A) 

-LN (0) 

v|/« (V) 


a < % > 

FF(%) 

v nr(™V) 

J« r (mA/cm 2 ) 

(1) 

10 

mm 

.6840 


14.440 

79.21 

484.32 

37.639 

15 

IB 

.7552 

.0000 

17.165 

81.70 

558.13 

37.644 


20 

1 

.7548 


17.288 

81.78 

651.58 

37.644 


25 

13 

.7497 

■ 

17.120 

80.96 

561.75 

37.644 

(2) 

10 

10 

.6831 

.0039 

16.587 

81.00 

544.77 

37.639 

15 

15 

.7534 

.0246 

17.157 

81.06 

562.30 

37.644 


20 

20 

.7431 

.1516 

14.117 

66.69 

562.35 

37.643 


25 

25 

.7285 

.2953 

8.755 

41.36 

562.37 

37.640 


Calculation was made for N A = 2 x 10 ,5 cm 3 , <t> w — 4.1 V, H — 250pm, 

e, = 1, N F =5x\0 n cm- 2 . 


Isc(mA) V nc (mV) EFF(S?) 



FF(55) 

100 

\ 90 

I 80 

\ 70 
\ 60 
I- 50 
\ 40 
30 

I 20 
: 10 
0 


Figure 1: The efficiency, the open circuit voltage, the short circuit current, and 
the fill factor for the MIS solar cell for different oxide thicknesses di 
with the new assumptions. 
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Isc(mA) V„ r (mV) EFF(%) 



FF(SS) 

100 

: 90 

: 80 

: 70 

: 60 

: 50 

: 40 

: 30 

: 20 

: 10 

1 0 


Figure 2: 


The efficiency, the open circuit voltage, the short circuit current, and 
the fill factor for the MIS solar cell for different oxide thicknesses d; 
with eq uili brium conditions and old tunneling factor. 


1(A) 



V(V) 


Figure 3: The I-V characteristics for MIS solar cells. 

1) dj = 1 5 A, equilibrium and old tunneling factor. 

2 ) dj = 1 5 A, non equilibrium and new tunneling factor. 

3) dj = 20A, equilibrium and old tunneling factor. 

4) dj = 20A, non equilibrium and new tunneling factor. 

5) dj = 25A, equilibrium and old tunneling factor. 

6) d, = 25 A, non equilibrium and new tunneling factor. 
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The effect of the third assumption is illustrated in Table 6. As the permittivity of the thin 
oxide changes from 3.9 to 1 while fixing the metal work <f> M , the hole tunneling current J pt 

increases; the electron tunneling current J nt decreases; the diffusion current J dn , the electron 
current through surface states J m , and the generation-recombination current in the depletion 
region increase. The increase of J dn , J ns , and causes a decrease in the fill factor. As 

the dark-current components increase, the efficiency, short-circuit current, and open-circuit 

voltage decrease. 

The effect of this assumption is relatively small when the equilibrium conditions are 
satisfied. When the surface-state densities or the doping concentrations of the substrate are 
changed to get nonequilibrium conditions, the effect of permittivity values according to the third 
assumption is noticeable (see Table 6). The nonequilibrium conditions mean a substantial 
decrease in surface potential vj/, which leads to a bigger increase in the voltage drop across the 
insulator. The dark-curTent components become bigger and the device is forced deeper into 
nonequilibrium when the permittivity e, changes from 3.9 to 1. The efficiency, the fill factor, and 
the open-circuit voltage will be reduced. The same effect of changing e, can be noticed with 
changing thin oxide thickness d t . For very thin thin oxide (d t <10A) or veiy thick thin oxide 
(dj > 30 A) the effect of the value of e, is significant because of the existence of the 
nonequilibrium conditions. But with 15A «/, <25A the near-equilibrium conditions occur and 
the effect of the value of e ; is negligible on the performance of the MIS solar cells. 
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Table 6. Effects of Different Values of the Oxide Permittivity 


e i 

para- 

meter 

Vs OO 

<1 >.(V) 

T1 (%) 

FF(%) 

V oc (mV) 

J sc (mA/cm 2 ) 

1 

+nr“>‘ 

♦m -4 - 7 

.7548 

.2798 

Mtid 

1 

m 

81.777 

60.624 

561.58 

185.21 

■ 

1 

3.9 

♦nr 4 - 1 

^nr" 47 

.7867 

.2425 

.0000 

.3282 

17.310 

2.880 

81.975 

55.485 

■ 

H 



1 Dj t =1.10 10 

Dj t =5.10 12 




N*=2.10 17 


N a =2.10 13 

N a =2.10 17 


7.318 

.537 


0 

.8043 

2 

.6232 


.6392 

.5914 




.0000 

.0001 


3.231 

2.649 


3.237 

7.022 



562.37 

42.79 


562.39 

468.30 


442.56 

521.77 



442.55 

38.229 

636.29 

32.476 



Nf=10 9 .660 .0000 16.745 81.219 547.77 37.637 

1 057 .0008 17.351 81.694 562.36 37.768 
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The effects of changing the values of parameters on the behavior of the MIS solar cells 
and MIS contacts are summarized as follows: 

1. The Effects of Changing Metal Work Functions <|> m 

The different values of work functions of metals, <f> m , affect the current-voltage 
characteristics of MIS contact under dark conditions. For <J> W less than a value of 
d> f<b <d> ) the semiconductor is strongly inverted. The dominant component of current is 

the diffusion current for reverse bias and for small forward bias. For higher forward applied 
voltage, the MIS contact is driven to depletion, and the current becomes larger. The main 
component of current is now due to the recombination in the depletion region. The hole and 
electron tunneling currents through the thin oxide are negligibly small. For <J> m > (j) mc , the MIS 
contact is no longer at strong inversion, except for very high reverse bias. The dominant 
component of dark current is now that passes through the surface state J ps . The diffusion 
current and the recombination current in depletion region become less important and can be 
neglected The tunneling through surface states dominates and controls the behavior of the MIS 
contact. The dark current J dark is approximately equal to the component that passes through the 
surface state J ps [3]. The I-V characteristics of a MIS contact for different 4> m are illustrated in 

Figure 4. 
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4. 4.3 «V 




For MIS solar cells the choice of the metal is essential. Figure 5 shows the effect of 
choosing metals with different work functions on the efficiency, the fill factor, the short- 
circuit current, and open-circuit voltage. For <j> m of less than a certain value (<t> mc ) we get higher 
efficiency, increasing <j> m decreases the efficiency. Our calculations for dark conditions [3] show 
that <i> = 4 3 V under zero terminal voltage corresponds to the onset of strong inversion for the 

parameters used in our model. For lower <t> w < <J> mc , the split between the electron and hole Fermi 
levels 4>, is approximately equal to qV, where V is the terminal voltage. When the MIS cell is 
under illumination, the incident photons generate electron-hole pairs. The electrons go to the 
surface. Because the surface is initially strongly inverted, the concentration of electrons at the 
surface is not affected too much by the light and near equilibrium conditions can be assumed. 


Isc(mA) V 0C (mV) 


EFF(%) 



FF(S5) 
00 

: 90 

! 80 

70 

60 

: 50 

\ 40 

\ 30 

20 

10 

«^0 


Figure 5: 


The efficiency, the open circuit voltage, the short circuit current, and 
the fill factor for the MIS solar cell for different metal work functions 

^m- 
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The current drawn from the cell can be written as J( — Jught Jdark 
Jdark = Jps +Jpt +J dn+ J rg Upt is the hole tunneUng current. J ps is the electron current 
through surface state. J dn is the diffusion current. J ^ is the generation recombination current 
in the depletion region.) For 4> m <4> mc , is approximately equal to the terminal voltage and 
j . is small. The main component of the dark current is J dn and the open-circuit conditions 
occur when J dn = J light . This leads to higher open-circuit voltage (V^), fill factor (FF) and 
efficiency (rt). For <j) m ><j) mc the tunneling current through surface states (J ps ) becomes the 
dominant current component. The open-circuit condition occurs when J ps = Ju g ht ■ ^ occurs at 
lower terminal voltage. The increase of 4> m increases J ps . As a consequence, and efficiency 
r| decrease and, hence, the fill factor decreases. ^ 

2. The Effect of Changing the Substrate Doping Concentration N a 

The doping concentrations of the silicon substrate affect the dark current of the MIS 
contact. There is a critical value N ac under which the I-V characteristic follows the ideal diode 
equation. This N ac is found to be around 10 17 c/n 3 . For N a < N ac the main components of 
dark current are the diffusion current J dn for reverse and small forward bias, and the 
recombination current in depletion region for reverse bias. The I-V characteristics of a MIS 
contact for different N a are illustrated in Figure 6. The main tunneling current is due to the 
minority carrier under the dark conditions. 
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The dependence of efficiency, open-circuit voltage, short-circuit current, and fill factor on 
the substrate doping N a is shown in Figure 7. The figure demonstrates that for a lightly doped 
substrate the efficiency is low. By increasing N a the efficiency, the open-circuit voltage, and fill 
factor increase until they reach the peak for the range N a = 10 15 to 10 17 cm" 3 Then, they start 
decreasing a gain For the short-circuit current, it decreases while the doping is increased because 
the lifetime of minority carriers decreases with increasing N a . Using a certain metal work 
function, the band bending % increases with increasing N a due to the increase of metal- 
semiconductor work function difference until a certain doping density ( N a —5x10 cvn ) is 
reached, which is called N ac . For N a <N ac the efficiency is increased due to the increase of 
band bending . In this case, the semiconductor surface is thermally inverted and equilibrium 
conditions occur. For N a > N ac the semiconductor charges at the surface decrease until the 
effect of surface-state charges prevails, which leads to a decrease in Tj until the inversion layer 
vanishes. For a similar reason like that for higher <j> m > <|> mc , the open-circuit voltage, efficiency, 
and fill factor become lower with increasing N a > N ac . 
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3. The Effect of Changing the Surface-State Density D; t 

The effects of the surface-states density D it depend mainly on the choice of other 
parameters (4> m , N a , (the thickness of the thin oxide), etc.). For choosing nearly optimum 
parameters, we get an inverted semiconductor surface. Increasing the surface-states density 
de cr ea ses the band bending and causes a decrease of the surface potential 'F f . As a result the 
dark current component J pt will increase and the minority-carrier tunneling current J nt will 

decrease. 

When D it exceeds a certain level which is called D itc , the effect of surface states as 
recombination centers disturbs the near-equilibrium conditions. The concentration of holes at the 
surface increases and the hole recombination current through the surface states ( J ps ) and the hole 
tunneling current ( J pt ) increase significantly. As a result, the dark current increases substantially 
with increasing D it . Therefore, the efficiency, the fill factor, and the open-circuit voltage 
decrease as shown in Figure 8. 



densities Dj t . 
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As D it increases further, the recombination current through surface states ( J ps ) becomes 
very large even at short-circuit conditions. For good quality semiconductors using hydrogen 

annealing, the surface-state density is in the range of 10 —10 states I cm ■ eV , and their effect 

on the performance of the MIS solar cells is reduced if the other parameters N a , N j- (the 
fixed oxide changes), dj ) are well chosen. 

4. The Effect of Changing the Thin Oxide Thickness dj 

In our models we have studied the effects of changing the oxide thickness dj which is an 
important parameter which affects the performance of MIS solar cells. Under dark conditions, the 
main components of current through the thin oxide in a MIS contact are J r g (for reverse bias and 
small forward bias) and d Dn (for large forward bias). The hole tunneling current J pi is high, but 
less than J I3W . As the thin oxide thickness d { increases, the tunneling current J pt drops 

substantially. This is shown in Figure 9. 

In the case of MIS solar cells, the effects of changing d { depends on the choice of the 
other parameters. Figure 1 shows the effect of dj . For thinner oxide the hole tunneling becomes 
easier, which means that the dark current becomes larger. The total current becomes smaller and 
the open-circuit voltage is reduced and the efficiency decreases. As the thickness of the oxide 
layer increases, the hole tunneling through the insulator decreases which leads to a decrease in 
dark current. So the efficiency and open-circuit voltage increase until a certain thickness of oxide 
d ic is reached. For dj>d ic , the electron tunneling becomes more difficult and higher 
concentration of photogenerated electrons exists in the semiconductor. As a result, an optically 
induced inversion layer is created and the split between electron and hole Fermi levels increases, 
so the dark current increases again. This causes a decrease in the open-circuit voltage and 
efficiency. When dj > d ic and dj reaches a certain value, the dark current becomes very large 
even at short-circuit conditions and the total current becomes very small (the short-circuit current 
J sc decreases rapidly and reaches almost zero) as shown in Figure 1. For more detailed 
discussions see references [3] and [4], 
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The most important parameters of the inverse layer MIS solar cells, for design 
consideration, include the doping concentration of the substrate and the thickness of the interfacial 
layer of the thin oxide in a MIS contact. Therefore, a lot of conclusions drawn from the study of 
MIS devices can be applied in the study of the inversion-layer MIS solar cells because MIS 
contacts are used in both cases. 

4.0 The Two-Dimensional Model for the Inversion-Layer Metal-Insulator- 
Semiconductor (MIS) Solar Cells 

A general purpose, two-dimensional, microelectronics device simulator using the finite 
element method has been developed [15]. This model is applied to the inversion-layer MIS solar 
cells. To model the MIS contacts the band to band tunneling currents and the current 
components due to the surface states need to be determined for the calculation of the dark 
current. The band to band currents have two components, one is the hole tunneling current and 
the other is the electron tunneling current. The hole tunneling current is 

Kfp-Evo & fm-Kyo 

Jp = " -<T ] W 

The electron tunneling current is 

Eco~Efin _ & co ~&fn 

J nt = 4T 2 Q(X„, d)[e—B^-e- * ] ( 4 ) 

where 0(x,<f,) = exp(-x 1/2 ^,), which is the tunneling factor 0(x,^,) = exp{f|A^<it:} [16] for 
J| Kfa=x m di, and the curves from Card [10,11] are used to obtain y}' 2 ^i A p “d A* are 
effective Richardson constants for holes and electrons, respectively. X«(X P ) is an effective 
potential barrier of oxide to electron (hole) tunneling into metal. 

The surface-states current components include the electron current through surface states 
and the hole current through surface states. The former can be expressed as 

= < 5 ) 
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The latter can be written as 


J ps ~ yDitVthO p\ftPs 0 ft )Pil 
In equations (5) and (6) 


E ti~ E W 

n x = N c e KT 

(7) 

Efr-Evo 

Pi = N v e ” 

(8) 

Kfh-Eco 

n s = N c e 

(9) 

E fp-Eyo 

p s = N v e 

(10) 

ft ~ T,+T, 

(11) 

n/5 n + PxO p 

ft0 ~ (n s +n l )c„+(p s +p l )o p 

(12) 



(13) 

(14) 

(15) 
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(16) 


J Dn ~ 


s£s. n 

L. P °F,(H) 


- 1 ] 


where 


F X {H') 


SnLn 

A, 


. ,H' ,H’ 

sinh k cosh — 

A» A 


(17) 


S„L, u H ' ■ U H> 

F 2 (H') = -2-72- cosh— +sinh— 

D„ A, Ai 

The generation recombination current in the depletion region can be written as 


(18) 


qrii 


Jrg 




2e 5 v|/, 




[e 


( W ,)/2 _ 


1] 


(19) 


All these current components are shown in Figure 10. In this diagram the Fermi levels for 
electrons and holes are split with energy difference equal to q$ s , the voltage drop across the 
insulator is A, the band bending is v, , the neutral level for surface states is E ss . 



Figure 10: Energy band diagram for MIS contact. 
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( 20 ) 


The voltage drop across the oxide can be written as 
a = E g +x a -b m -r p -v s -v 

where V is the terminal forward bias. 

The dark current in this simulation is obtained by 

Jdark - Jfit + Jpt ~ J ns + ^ ps 
where J nt is the electron tunneling current, J pt is the hole tunneling current, is the electron 
current through surface states, and J ps is the hole current through surface states. 

The effect of light is treated as a generation term which is given by 

g(x) = J* [1 ■ - /?(X)]a(X) N mfr) exp(-a (\)x)dk (22) 

where N m (X) is the mass (m) photon flux spectrum, a(X) is the silicon absorptivity, and R(k) 

is the optical reflection of the cell's front surface at the wavelength (X) [18], 

Figure 11 is a partial mesh of the inversion layer MIS solar cell. The equations for 



Figure 11: Partial mesh of the inversion-layer MIS solar cell. 


24 



tunneling currents and for surface state currents (Eqs. (3) to (6) ) are applied to the thin oxide 
area which is between the metal and the silicon in the cells. The diffusive and drift current 
densities are calculated by Two-Dimensional General Purpose Finite Element Microelectronic 

Device Simulator [15]. 

The simulation of the inversion layer MIS solar cell using this two-dimensional 
microelectronic device simulator is a very difficult task. The very thin oxide layer between metal 
and semiconductor is about 20 A, which is much smaller than the cell dimension used in this 
simulation ( 300 \xm x 80 \mt ). This is very difficult for the numerical calculation. 

The values of the parameters used in this simulation are shown in Table 7. The simulation 
results of the short circuit currents for different thickness of the thin oxide are shown in Table 8. 
The simulation results do not agree with the experimental values [19]. The problems come from 
the simulation of tunneling currents. More time is needed to have a further study on the effects of 
the inversion layer, the surface states, and the tunneling in the MIS contacts. 


Table 7. Short Circuit Currents for Different Thickness of Tunneling Area [15] 

(X is the barrier potential for tunneling area, 8 is the thickness of tunneling area) 


tunneling area thickness 

X m 8 (eV 1/2 A) 


15 A 

2.7 

no convergence 

20 A 

7.5 

0.337 

30 A 

12.5 

0.124 
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Table 8. Parameters Used in Simulation 


PARAMETER 

VALUE 

UNIT 

D;t 

6.0 Ell 

states/cm 2 /eV 

Vth 

0.0256 

V 

Na 

2.0 E15 

/cm 3 

Nf 

5.0 Ell 

cm* 2 

4>m 

4.1 

V 

di 

20 

A 

An* 

112 

A/cm 2 .K 2 

An 

32 

A/cm 2 K 2 



E^ 

0.33 

eV 

n; 

1.45 E10 

/cm 3 

o n 

1.0 E-15 

cm 2 

2d 

1.0 E-15 

cm 2 


5.0 Results and Discussion 

Many 1 cm 2 inversion layer MIS solar cells have been fabricated. Most of them have low 
efficiencies, small short circuit currents, and small values of filled factors. There are eighteen 
1 cm 2 H/MIS cells which have the efficiencies between 6% and 7%, and there are four 1 cm 2 
cells which have the efficiencies between 8.06% and 8.9%. The highest short circuit current is 
30.400 mA while the largest open-circuit voltage obtained is 0.575V. 

The performance of the 19 cm 2 cells is poorer than that of the 1 cm 2 cells. The best 
19cm 2 IL/MIS cell has = 0.536 V, I sc = Sl\mA , FF = 0.394, and 5.374% efficiency. 

The growth of the thin oxide in the IL/MIS solar cells has been successful. The thin oxide 
in the range from 18 A to 25 A was successfully grown. Our analytical model for Metal- 
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Insulator-Semiconductor (MIS) devices is very successful to predict the effect of the thickness of 
the thin oxide on the efficiencies of the MIS solar cells, which is applicable to the IL/MIS solar 

cells. 

A lot of efforts have been made to develop a two-dimensional model for the IL/MIS solar 
cells. Unfortunately it has not been successful. More work is necessary to develop the two- 
dimensional simulator for these photovoltaic devices. 
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